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ABSTRACT 

In developing turbulence models, different authors have proposed various model con- 
straints in an attempt to make the model equations more general (or universal). The most 
recent of these are the realizability principle (Lumley 1978, Schumann 1977), the linearity 
principle (Pope 1983), the rapid distortion theory (Reynolds 1987) and the material indif- 
ference principle (Speziale 1983). In this paper we will discuss several issues concerning 
these principles and will pay special attention to the realizability principle raised by Lum- 
ley (1978). Realizability (defined as the requirement of non-negative energy and Schwarz’ 
inequality between any fluctuating quantities) is the basic physical and mathematical prin- 
ciple that any modeled equation should obey. Hence, it is the most universal, important 
and also the minimal requirement for a model equation to prevent it from producing un- 
physical results. In this paper we will describe in detail the principle of realizability, derive 
the realizability conditions for various turbulence models, and propose the model forms 
for the pressure correlation terms in the second moment equations. Detailed comparisons 
of various turbulence models (Launder et al. 1975, Craft et al. 1989, Zeman and Lumley 
1976, Shih and Lumley 1985 and one proposed here) with experiments and direct numeri- 
cal simulations will be presented. As a special case of turbulence, we will also discuss the 
two-dimensional two-component turbulence modeling. 
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1. Introduction 


It is commonly accepted that the Navier-Stokes is the basic equation governing tur- 
bulence. Hence, without any model approximations, direct numerical simulation (DNS) 
should be an ideal approach to study all the aspects of turbulence. For example, using DNS 
we may study the dynamics of various turbulence structures; the mechanisms of energy 
transfer; and the interaction between different size of eddies; and we may also look at the 
various one-point and two-point statistics. Unfortunately, DNS is limited to relatively low 
Reynolds number and simpler flows due to a wide range of length and time scales that are 
present in a turbulent flow and due to today’s computational limitations (both computer 
hardware and numerical methods). Therefore, at present, DNS is far from computing 
the practical engineering problems. For solving general engineering problems, turbulence 
modeling must be used in conjunction with the CFD. The history and the literature of 
turbulence modeling is described in detail by Lumley (1978,1983), Reynolds (1976), Laun- 
der (1989) and Rodi (1988). All the model equations were designed to approximate the 
Navier-Stokes equations. The model coefficients (or constants) have been tuned against 
some basic turbulent flows. The accuracy of any model equation (if it works) varies from 
flow to flow. In addition, we must use them with caution because any model equation will 
miss some of the physics or violate some of the properties of the Navier-Stokes equations 
because of the approximations made. Results in such cases may be misleading in studying 
unknown turbulence phenomena. More severely, some of the model equations may end up 
with unphysical results, such as negative turbulence kinetic energy, etc. 

For constructing turbulence models, various model constraints have been proposed 
by different authors in an attempt to make the model equations as general (or universal) 
as possible. Most recent of these are, the realizability principle (Schumann 1977, Lumley 
1978), the linearity principle (Pope 1983), the rapid distortion theory (Reynolds 1987) and 
the material indifference principle (Speziale 1983). However, some of the above mentioned 
principles are not universal. For example, the material indifference principle is not valid for 
general turbulence in which the fluctuating velocities are three dimensional. The material 
indifference is valid only for strict 2D-2C (two dimensional and two velocity component) 
turbulence which rarely occurs or probably does not exist in nature which is described by 
three dimensional space. Speziale et al. (1987), in their DNS experiments, found that rapid 
rotation had no influence on the isotropy of the 3D-3C homogeneous isotropic turbulence. 
So far, there is no evidence which shows a 3D-3C fluctuating velocity field will evolve to 
a strict 2D-2C state under any conditions. Therefore, the material indifference (or frame 
indifference) principle need not be considered in 3D-3C turbulence modeling. In addition, 
as we will see later, for strictly 2D-2C turbulence the frame indifference principle principle 
will be automatically satisfied and does not provide any additional or new constraints 
for turbulence models. In another study, Reynolds (1987) also showed that in 3D-3C 
turbulence modeling, at the 2D-2C limit, the frame indifference principle does not provide 
any new model constraints and it also seems to be irrelevant to the modeling of the effect 
of the rotation on turbulence (see the RDT described later). 
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The principle of linearity is also not universal, as it holds only for passive scalars. 
It was pointed out by Pope (1983) that since the governing equations for passive scalars 
are linear, the modeled equations for these quantities should possess the same property. 
However, some of the models for the pressure correlations are nonlinear. The calculations 
of slightly heated turbulent jets (in which the temperature is a passive scalar) with such 
models (Shih et al 1990) show that the violation of this principle does not cause problems 
since the results agreed well with the experimental data (also see Lumley 1983). 

Rapid distortion theory (RDT) is a convenient tool for studying the asymptotic be- 
havior of turbulence under the rapid distortion of the mean flow. For example, effect of 
the mean flow (say strain rate) on turbulence. For some simple homogeneous flows RDT 
can provide an analytical solution (Reynolds 1987), hence enabling us to calculate the 
model terms of interest in the second moment equations. One may naturally ask if we 
can use these results to calibrate the model coefficients appearing in the corresponding 
model terms. Unfortunately, the answer is not very positive. RDT equations are valid 
only for very high mean strain rate or very low turbulent intensity and in both these cases 
the turbulence is in a very non-equilibrium state with the mean flow field. On the other 
hand the second order modeling is based on the assumption that the turbulence is more 
or less in an equilibrium state and is mainly characterized by one velocity and one length 
scale (Lumley 1983). Therefore, RDT results are probably not appropriate for tuning the 
coefficients of second order models. An interesting illustration of this point is made by 
Shih et al. (1990) who showed that a model constant in the linear form of rapid pressure 
strain model (Launder et al 1975) determined by matching RDT results is much worse 
than the one determined from measurements. However this does not mean that RDT is 
useless in turbulence modeling. On the contrary, RDT can provide the asymptotic behav- 
ior of turbulence at the RDT limit. In many cases, we may obtain analytical expressions 
which explicitly show the relations between turbulence and the mean flow field. For ex- 
ample, the RDT results (Reynolds 1989) for pure mean rotation show a clear effect of the 
mean rotation on the anisotropy of turbulence (bij = u x ujj q 2 — 3) through its second 

invariant. However, all the existing rapid pressure strain models (including the one which 
satisfies the material indifference at the 2D-2C limit) fail to exhibit any effect of rotation 
on the second invariant of 6,y. In this case, RDT will certainly guide us in removing this 
deficiency from the existing models. 

Realizability (defined as the requirement of non-negative energy and Schwarz’ inequal- 
ity between any fluctuating quantities) is the basic physical and mathematical principle 
that any governing equation should obey. Hence, among all the above mentioned model 
constraints, realizability is the most universal, important and also the minimal requirement 
for the turbulence model equations to prevent them from producing unphysical results. Re- 
alizability can be applied to various turbulence quantities, for example, Reynolds stresses 
ttjttj, scalar fluxes u,6 and the triple correlations. For one-point second moment equations, 
realizability for TTpuJ and urf is the most important. In this paper (section 3, also see Shih 
and Lumley 1985) we will describe the principle of realizability, derive the realizability 
conditions for various model terms, and propose the model forms for pressure correlation 
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terms in the second moment equations. The comparisons of various turbulence models 
with DNS and measurements will be presented in section 4. However, first we will discuss 
the 2D-2C turbulence. 


2. Two-Dimensional Two-Component Turbulence 


Reynolds and Grau (1988) argued that if one component of velocity is suppressed and 
turbulence becomes 2C, the turbulence will not stay in the 2C state unless the turbulence 
also becomes two dimensional (2D). In the DNS experiment, Speziale et al (1987) were 
expecting that the homogeneous isotropic turbulence would become anisotropic under the 
action of rapid rotation. However, their results did not show any such trend. Can a 
wall bounded inhomogeneous turbulence become strictly 2D-2C turbulence under rapid 
rotation? In a three dimensional space, it seems unlikely that the instantaneous turbulent 
velocity field will become strictly or exactly 2D-2C turbulence. However, if we confine 
a flow to a two dimensional space, as in the “thin soap sheet film”, then it may be possible 
to obtain a 2D-2C turbulence. Numerically it is always possible to simulate such a special 
“turbulence” . 

In this section, we will derive a 2D-2C turbulence model for the rapid pressure corre- 
lation term in the Reynolds stress equations (also see Reynolds 1987) and will show that 
the frame indifference principle does not provide any model constraints. 

The Reynolds stress equation can be written as 


— UiU q = ... + T iq - 2vu i>k u qyk (2.1) 

where 

% = 2 U rJ [M ijn + M, ipi ) (2.2) 

r 

corresponds to the rapid part of pressure correlation term, and q 2 = uiul. For homogeneous 
turbulence, Mij pq can be expressed as 

Mii " = 7 LQ E “ ik)d * k (2 ' 3) 


where k{ is the wave number vector and Ejj is the spectrum tensor. 

As usual, we assume that Mij pq is a function of the anisotropic tensor bij, and then 
the most general tensor form of Mij pq (which satisfies the symmetry in indices) will be 
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(2.4) 


where, 


— Q'l&ij&pq “I” ^2 (^ip&jq “I" ^iq^jp') d" "f” ^4 b pq bjj 

~f“ Q’b^iph jq “H ^iqbjp ~b bjpb{ q d" ^jq^ip') 

+ ^6 bijb p q + aq(bi p bjg + bi q bj p ) 


Sii = 2 

II = —bijbij/2 


(2.5) 


Because of the Cayley-Hamilton theorem in a 2D space, b?- + IIS l} = 0, the higher powers 
of bij are not required. Seven coefficients a, are, in general, function of the invariants of 
bij. From the continuity and definition of we have 


M pjpq = 0 

Mijpp — bjj -f- $ijf 2 

These give us the following relations between the coefficients: 


( 2 . 6 ) 


a\ + 3<22 — J/ag — Ilctj = 0 
a 3 + a 4 + 4a5 = 0 
2,q>i “I - 2 g 2 — 2Z/0E7 = 1/2 
2a 4 “I - 4as — 1 


The coefficients a 4 , a 2 , 03 and a 4 can be solved in terms of 05 , ag, and oj: 

a x = 3/8 + II(a 7 — «6/2) 

<22 = — 1/8 + IIa 6 /2 
<23 = — 1/2 — 2ag 
<2 4 = 1/2 — 2(25 

Therefore, the tensor Mij pq becomes 

Mijpq ™ g(3(5|ji6pg ^ip^jg ^ig^jp) 4~ 2 (^pg^ij ^tj^pg) 

+ ^6 [^ij bpq ~ Il&ij&pqfe + II(8{p8jq + 8iq8jp)/ 2] 

4“ &7 (^ipbjq “j - ^ig^jp 4“ -^T^ij^pg) 
where as drops out because of the following zero tensor: 

8 ip ^jq 4“ 8 i q b j p 4" 8jpb{ q 4* ^jg^ip 2<$,’j b pq 2 8 pq bij = 0 


(2.8) 


(2.9) 


(2.10) 
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( 2 . 11 ) 


Using (2.9), the Eq.(2.2) becomes 

T 

-==■ = Siq + 4(a 6 + 2a-j^bijbpqS p j -f- 2 Cl q jbij + 2 flijbjq 

r 

where we have used the following relation: 

4// Siq + [4 biqbjp — 2{bijb p q + bi p bjq)]Spj = 0 (2.12) 

Equations (2.10) and (2.12) can be easily checked in the principal axes. 

For the 2D-2C turbulence, the turbulence equation should be frame indifferent 
(Speziale 1983). If we write Eq.(2.1) in a non-inertial coordinate system (designated by *) 
which is rotating with an angular velocity Ojt, we obtain 

U ' U q = ^(Up,j ) \Mijpq ~b Mqjpi] q 2, UUi,kUq t k • 

2& jxkQk'U’ jHq 2c jqk& k^ 1 j'U'i 

£pj k \Mijpq "b ^^qjpi\ 9^ 

In order for this equation to be frame-indifferent, the last two lines must cancel each other. 
This gives us the necessary constraint on the model form of M ljvq and the result is 

04 - a 3 = 1 (2-14) 

However, this condition is already satisfied by the coefficients (2.8) and hence the frame 
indifference principle does not provide any additional constraints. Therefore, in general 
there is still one undetermined coefficient a^-\-2a-j in (2.11) and it must be determined from 
experiment, direct numerical simulation or some other source. However, we may explore 
its limiting value by using 1C condition. The condition of Mij vq at the limit of u\ = 0 
gives M\j pq = 0 for any j, p and q (because Eij — 0). This condition gives 

a 6 + 2a7 = — 1 if bn = — 1/2, (or uj = 0) (2.15) 

However, we must note that in general 2D-2C flows, the value of «g + 2a 7 is unknown and 
not necessarily —1. The limiting value of the above coefficient can also be obtained by 
using the realizability conditions (3.4 a,b,c) described in the next section. Interestingly 
enough such an analysis also gives a$ + 2a- = — 1. 


(2.13) 


3. Three-Dimensional Three-Component Turbulence. 


It is commonly thought that the turbulence in a three dimensional space is always 
three dimensional (Tennekes and Lumley 1972) with three components of velocity. How- 
ever, there are situations where the turbulence is considerably suppressed in a certain 
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direction. For example when turbulence approaches a solid wall, the fluctuating velocity 
component normal to the wall will be mostly suppressed and the turbulence will become 
3D-2C. However, during this process, all the turbulent energy components should remain 
non-negative and Schwarz’ inequality should always be satisfied. This is the concept of 
realizability. In this section, we will discuss realizability in detail and discuss the model 
development of pressure correlation terms in the second moment equations. 


3.1 Realizability Concept. 

In the second order closure, model equations for the turbulent stresses ujuj and scalar 
fluxes 6ui are proposed. Realizability for both of these quantities has been discussed 
in Lumley (1978). To see the realizability concept, let us consider any two fluctuating 
variables a and b, and their correlation matrix: 

(a 2 ab ) 

( ab b 2 J 


The simple fact that a 2 > 0, b 2 > 0, and that Schwarz inequality a 2 b 2 — ab > 0 should 
be satisfied by any set of equations governing the evolution of a 2 , b 2 , and ab is called 
realizability. A set of equations violating this principle could produce unphysical results. 
We will show that these three realizability constraints translate to the single requirement 
that the eigenvalues of the correlation matrix should remain non-negative. We will now 
show that if any one of the above three in-equalities becomes an equality, then one of 
the eigenvalues of the matrix will vanish. Therefore, by ensuring these eigenvalues from 
becoming negative will automatically ensure a 2 , b 2 , ( a 2 b 2 — ab ) from becoming negative. 
This implies that while constructing governing equations for a, b and their correlations, 
we need to check only the eigenvalues of the correlation matrix. The non-negativity of the 
eigenvalues of the matrix can be easily seen from their expression: 

Ai, 2 = (1/2) {(^ + ^) ± [& + ¥)* - 4(^2 ft 2 -^ 2 )] 1/2 } (3.1) 

Eq. (3.1) shows that A 2 will vanish if one of the following conditions occurs: 

~a 2 — 0, 
b 2 =0, 
b 2 b 2 - ~ab 2 = 0. 

Now, let us take the derivative of (3.1) with respect to time t, 

(Ai j2 ),t = (l/ 2 M( a2 ,t + b 2 >t ) 

± ±b\t) ~ 2[(^), t 6 ? + ^~ (¥), t - 2ab (ab),,] (3-2) 

[£? + 62) 2 -4(0* b 2 -ab 2 )}'/ 2 
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To preserve the positive A2, we must require 


(A 2 ),t >0, if A 2 — 0 

(3.3a) 

or 

fa 2 b 2 — ab 2 ] < 

i — = — = r-hi >0, if A 2 = 0 

a 2 + b 2 

From this necessary condition, it is straight forward to show that 

(3.36) 

a 2 }t >0, if a 2 = 0 

(3.4a) 

P,t >0, if ft 2 = 0 

(3.46) 

( a 2 b 2 — ab 2 ) it > 0, if a 2 b 2 — ab 2 — 0 

(3.4c) 


It is clear that (3.4) is the necessary condition for the quantities a 2 , b 2 and a 2 b 2 — 
2 

ab to remain non-negative, therefore, the condition (3.3a) is the necessary condition for 
realizability. To find the sufficient condition for realizability, we must examine the higher 
order derivative of the eigenvalue when it vanishes. For example, we may require 

(A 2 ),„>0, if A 2 = 0. (3.5) 

For more details on this aspect see the review paper by Lumley (1983). 

Similar arguments can be applied to the turbulent stress tensor ujuj and to the cor- 
relation matrix formed with a scalar 9 and velocities Ui. Realizability for tensor TTiuj 
is straight forward, that is, its eigenvalues must be positive, which will ensure positive 
turbulent energy components and the Schwarz inequality between fluctuating velocities. 
Realizability for the correlation matrix of scalar and velocity is described by Eq. (3.1) and 
(3.3) and is not as simple as for the Reynolds stress. Lumley has suggested a more general 
quantity Dij consisting of turbulent stress and flux: Dij = 0 2 TTJuj — 6ui 9uj. It can be 
shown that realizability for the correlation matrix between the velocity and the scalar will 
be ensured provided that the eigenvalues of Dij are positive. 

It is interesting to note that if u Q — ► 0 then the a direction will become one of the 
principal axes of uju], and u a u a (no summation on Greek indices) becomes one of the 
eigenvalues of UiuJ. At the same time, u a — > 0 will make D aa — » 0, and the a direction 
also becomes one of the principal axes of D tJ . D aa is then a vanishing eigenvalue of D tJ . 
However D aa will also vanish without u Q — ► 0 if the perfect correlation between / and 
u a is reached. Therefore, we conclude that if one of the eigenvalues of TtiuJ, say u a u a • 
vanishes, then D aQ will become a vanishing eigenvalue of the tensor Dij. However, the 
reverse is not true i.e. if D aa — > 0, it does not mean that u a u a (or 6 2 ) must vanish. This 
is important in the modeling of scalar turbulence, as it indicates that x’ealizability based 
on TTiuj alone is not enough and we must also consider the joint realizability based on Dij 
to ensure that the model equations do not produce unphysical results. 
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3.2 Realizability Conditions for Second Order Closure Models 

The ideas discussed above can now be applied to the second moment equations to 
obtain the realizability conditions. The general form of the model equations for the second 
moments can be written as follows 


(3.6) 


(3.7) 

(3.8) 

where we have followed Lumley (1978) for the pressure transport model and in combining 
the anisotropic part of the dissipation terms with the slow part of the pressure correlations. 
The molecular diffusion terms in the above equations have been neglected. The definitions 
of the rapid pressure correlation terms X p j q i and X,jk] the buoyancy terms Ykij and Yj*; 
and the slow terms $,j and are listed in appendix A. In the above equations is the 
buoyancy vector, 0 is the mean scalar and 9 is the fluctuating scalar. 


— = - [ UiUjUk - C{q 2 UiSj k + q 2 UjS ik ) ],* 

Uj,k "t" UjUk U{ f k) d~ 2Up iq (Xpjqi A piqj') 

- PkiXkij + Ykji) ~ (*y + 2Sij/3)e 

+ Piduj + f3j9ui 

= - [6uiu k - C'9q 2 6ik],k - (Ou k U iik + «iW0,fc) + M 2 
+ 2 U ith X ii k-foYik-*i*f? 

—9 2 = — ( 9 2 Uk),k ~ 2 9uk 0,jt - 2e# 


To impose realizability into the above equations we may write the equations of uJvTj 
and Dij in their principal axes, i.e., the equations for u a u a and D aa . Then in order 
to ensure that the eigenvalues of u^Ua and D aa remain positive during the evolution of 
turbulence, the following constraints must be satisfied 


D 

Dt 

D 


D 2 


D 


Dt 2 



->• 0 

if U a U a ■ 

) 

t Daa ' 

0 

if D aa 

^U a U a 

IV 

c 

if U Q U a 

2 



To ^ act 

> o 

if D CIO 


(3.9) 


(3.10) 


0 


Substituting (3.6)-(3.8) into (3.9) we obtain the following necessary realizability conditions 
for X p j q i, Xijk , Y{kj and Y, k ■ 


Up t qXp QqQ 


if u 


-> 0 

(in p.a. of Rij ) 
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-» 0 


U Piq (9 2 X pQqQ - 9u a X apq ) 

if D aa -*• 0 (u 2 a ± 0) (in p.a. of Dij ) 

U Ptq 9u a X a Pq * 0 

if D Qa ,u 2 Q -* 0 (in p.a. of Dij,Rij) 

Y n aa 9u a ► 0 

ifD aa ,u 2 a -> 0 (in p.a. of D tJ , R { j) 


Yy aa ( 7 / a ) 0 

if D aa ,u 2 -+ 0 (in p.a. of Dij,Rij ) 


* 2 n 

cm - 9u a Y ka (any k ) — ► 0 

if D aa -* 0 (u 2 ^ 0) (in p.a. of Dij) 

0u a Y aQ — 9 2 9u a — > 0 

if D aa ,u 2 a -+ 0 (in p.a. o/ Dij,Rij) 

9 u q Y^ q (7 ^ a) -*0 
if D aa ,ul-> 0 (in p.a. of Dij,Rij) 


($ af _+2/3)e ->F a 
if a 2 — » 0 (in p.a. 0 / i?,j) 


29u a $ a e/q 2 - 9 2 ($ aa + 2/3)e - 2e*«7 - F] 

if Daa -* 0 (a 2 ^ 0) (in p.a. 0 / D.j) 

2 Ih^Qat/q* - 2e e uT + ~9 2 F a -» F] 

if D aa , a 2 -> 0 (in p.a. of D^, R tj ) 

(u 2 a u k -2Cq 2 u a S ak ), k -» 0 

if u 2 -» 0 (in p.a. 0 / i?,j) 

2 9u a (9u a u k - C'q 2 9 6 ak ) tk 
- 9 2 (u 2 u k - 2Cq 2 u a 6 ak ), k -u 2 a (9u k ) ik -» 0 
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if D aa 


0 (ul ± 0) (in p.a. of Dij ) 


29u a (9u a u k - C'q 2 9 8 ak ) t k - ul(6u k ) >k -> 0 
if 0 (in p.o. of 


These constraints will be used in the next section in constructing the models for 
pressure- correlat ions . 


3.3 Modeling of the Rapid Pressure-Correlation Terms 

First let us look at the rapid terms associated with the mean strain. The symmetry 
property of these terms requires 


Xpjqi — Xjpqiy 


Xpjqi — Xpjiq, X{j k — Ji-jilc 


If we assume that X p j q i is a function of uiu] only and that X x j k is a function of both 
uju] and 9u,, but linear in 9u{ , then the most general forms of X VJqi and X{j k satisfying 
the above symmetry constraints are 

Xpjqi / f 2 = Otl&qi^pj 4" &2(bpq8ij 4" bqj8pf) 4" ^38 qibpj 4" ^4^pj^qi 

4" Qb(8 pqb{j "f 8ijb pq 4“ 8j q b px 4" ^pibqj') 4" ^6 bqib p j b ^7 ^pjbqi 
+ a 8 (^pq^ij "b ^ijbpq 4" 8 jqb px 4" ^pt^qj) 4" OlQ^qi^pj 
*b ^lo(bpqb{j + bqjbpf) 4* Otllbqibpj 4~ 

+ <*i3(b P qb?j + bijb 2 pq -b b q jb 2 pi + b p {b 2 -) 

+ a 14 bqibpj -b <*ib(b 2 pqb 2 j -b b 2 jb pi ) (3.12) 


X ik j = Pi8i k 9uj ~b $ 2(^1 j^'Uk + 8j k 9uj) -|- /3 3 bi k 9uj 

-b (3i(bij9u k 4” bj k 9u{) 4" fib(8ijb k p -b 8 k jb{p')9ii p 

4" $q 8 i k b j p9u p -j - /^qb{ k bjp9up 4~ @&(bijb k p 4~ b k jbjp')9up 

4- ^b 2 xk 9uj 4- /?io (b 2 j9u k 4- b 2 k 9ui ) 4- /?n 8i k b 2 p 9u p 

4- /?i 2 (^ij b k p 4* 8j k b ip )9u p -(■ fti3b{ k b j p 9u p 4* Pi4bi k bj p 9up 

+ Pib(b 2 jb kp + b 2 k jbip)9u p + /?i6 (bijb 2 kp 4- b kj b 2 p )9u p 

+ Pnb 2 k b 2 p 9u p 4- P\t,(b 2 jb 2 kp 4- b\-b 2 ip )9up (3.13) 

Where a’s are, in general, functions of the second and third invariants (i.e. II and III) 
of the anisotropic tensor 6,j and /?’ s are functions of the invariants formed by b XJ and 6u t 
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in addition to II and III. From the definition of the rapid terms and incompressibility, 
following properties are obtained: 


Xppqi — UqU j, Xpkqk — O 5 X Hk — ^ikk — 0 (3.14) 


Using these properties and the realizability conditions given in the last section we can 
determine the limiting values of all the model coefficients provided we ignore the terms 
with the fourth powers of b{j (which is a reasonable assumption since the magnitude of 
bij is always less than unity). The final expressions of the X p j q i and Xijk reduce to the 
surprisingly simple forms and are given by 


Xpjqi /(p fipqfiij fiqjfopi) ^i^qi^pj j ^q i ) 

11 4 

4“ Q>\[£pqbij “h Sijbpq “(" &qjbpi Sptbqj ~^~^qi^pj 
0*2 i^^pjbqi Qbpqbij Qbqjbpi + bq{bpj^j 


X{jk — ^bijOlLk 
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(SifcBlLj -f~ ^jk^^i) 


+ C D\ bij @Uk 

C E)2{pik$u j “I - bjk&y>i ) “l - C £)^S ijbki^y^ i 
where the limiting values of the coefficients are 

= “To’ “ 2 = To’ Cdi = To’ c “ = “To’ Cm = 5 


(3.15) 


(3.16) 


The last line in both (3.15) and (3.16) represents the non-linear contribution, and 
if neglected, the linear models used by various other workers will be recovered. It is 
important to note that the above values of the coefficients ai, a 2 , Cdi, C D2 , and Cm 
are their limiting values at the realizability limit, i.e. when u 0 n a , D aa — ► 0. For general 
turbulent flows and D aa are not zero and hence the values of the coefficients may 

deviate from their limiting values. As was pointed out before, in general a’s are functions 
of II and III, and Cd ' s are functions of the invariants formed by the anisotropic tensor 
and the heat flux. Some guidance can be obtained by inspecting the following two useful 
parameters (see Lumley 1978, and Shih and Lumley 1986) 

F = 1 + 27/// + 97/ (3.18) 

ft = 9 4- f 4 + f (3.19) 

where 

u = -\bl, m = ±bl 

6 2 TTpuJ — dui duj 
u 6 2 — dui dui 
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It can be shown that both F and Fd are bounded between 0 and 1, and particularly 


F — ► 0 when u a u a — ► 0 
Fd — * 0 when D aa — ► 0 

By using this information, it is convenient to write 

a ‘ = ~To (1 + AFa) (3 ' 20) 

»2 = ^(1 + BF°) (3.21) 

Cdi = ^ + Cl Fj (3.22) 

Cm = ^ + C, fy (3.23) 

Cm = 7 + C,F; (3.24) 

where A, £?, C \ , C2 and C3 are adjustable constants but a is not as arbitrary as it might 
seem at first look. The condition (3.10) suggests a = ^ (see Lumley 1983, or Shih and 
Lumley 1986). In the limiting case these coefficients reach their limiting values. Shih et 
al (1985, 1987) took A — 0.8 and B = 0.0 which fit the DNS data quite well. However, 
it seems more judicious to use homogeneous shear flow experiment to obtain the values 
of these coefficients (this is not done in the present study). Shih et al (1990) set C\ = 
C2 = C3 = 0 in their computations. The present study indicates that the full form of 
Cdi,Cd2 and Cd 3 should be used. The experiment of Tavalouris and Corrsin (1981) 
suggests Ci = 1.8, C2 = —1.8, and C3 = 4.5. 

Now returning to the modeling of the buoyancy terms Yjt,j and Y, k , we assume that 
these are functions of the anisotropic tensor and the heat flux. Their minimum forms 
required to satisfy realizability are as follows 

Yikj = 0\8{ k 9uj + 0 2 (8ij9u k + 8j k 9u{) + 03 b{ k 9uj 

-f- 0\{bij9u k -J- bjkdui) ftsihijbfcp -(- hkjbip')9up 

"t 0< 6btkbjp9up -f- 0-jbikbjp6up ~b 0&{bijb) C p A b k jb{ p ^9ii p 

+ 09b 2 k 9uj + 0io(b?j9uk + b 2 j k 9ui) 4 - 0n8i k b 2 p 9up 

+ 0\2{8ijb\ v + 8j k b 2 ip )9up (3.25) 

Yi k = 7i 9' 2 8{ k T i 2 e 2 b ik b 73 9u{ 9u k 4 74 (9u k b{p9iip 4 9u{h k p9v , p ) 

4 + 76 b ip b kg 9u p 9u q 

4 77 {8} p 9up 9 u k + b\ p 9up 9ui ) (3.26) 

The above forms already satisfy the following symmetry conditions 

Y klj = Y ikj , Y tk = Y kt 
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By definition, and incompressibility, these tensors have the following additional properties: 

Y kkj = e^, Y kk =¥ 

Using these and the realizability conditions given in the last section, all except one coeffi- 
cient can be determined exactly. Their final forms are 

ft =|(i - ft) 


ft : 2 — 

A = 


(A - 1) 

3 

(i-nA) 


A =A - 1 


A 

A 


- (127/ + 7)rj 


6II9 2 q 2 — 10r 2 — 3 6/ Jrfrr 

(2 — 4 A) 

3 


6JJ + ^3JJ ^ 

A =3 A 

ft =67/-«377> 

Ao =o 
Ai = — 2A 
A 2 =o 

(6r 2 + 1O0 2 g 2 ) 


7i 


72 


90 2 g 2 

[18(rj + ¥ q*)rbr + rj(7F g - 15rj) + (36JJ - lO)^ ^) 2 ] 

90 2 g 2 (r£ — # 2 g 2 ) 

(g 2 /9)(3 rbr + 6/Jr 2 — 2II9 2 g 2 ) 

~~ (r£ - tf 2 

ri(6/J + 28) - fl 2 ?(6JJ + 7) 7(rj) 2 

18//(t£ — ¥ g 2 ) (6/J)? 2 ¥(rl-¥ q 2 ) 


-A 

+ 76 


rj(12/J + 20) + 108/Jrk + fl 2 g 2 (108// 2 - 21// - 5) 
9 II(r ±-¥ ?) 

5(3// + l)(r!) 2 
3/J0 2 ¥) 

r r 2 q 2 (18JJ + 4) + 9q 2 rbr — 8 2 g 2 (6// + l)g 2 
+ 76 ^ 9(ri - ^ g 2 ) 
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4 —] 

/lO o \ J 


30 2 (ri-^ 2 q 2 ) 

—2 12 r6r — 5r 2 + (127/ — l)fl 2 g 2 

73 ==- + P5 


-76 


q 2 (r±-9 2 q 2 ) 

6/Jd + 3r6r - 2W T 2 


3(li — 9 2 q 2 ) 


74 =-= 


-3ft 

2 


7 _ (18/7 — 5) ,g 2 r 2 . 

75 “6/7 + 3/7 + 76 

76 = undetermined coef ficient 

77 = 0 


where = 0u, 9u{ , r6r = 6 tJ - Quj. These expressions may seem a bit complicated 
but they ensure joint realizability. Since 76 is undetermined it was taken as zero in the 
calculations to be presented later in this paper. As an illustration, Appendix B lists the 
values of all the above coefficients as calculated from the buoyant plume experiment (1990). 


3.4 Slow Terms 

Lumley (1978) has suggested the following expression for the slow term of the pressure- 
strain correlation 

*ij = fibij + 7 (bikbkj + 2IISij/3) (3.27) 

where f3 and 7 are, in general, functions of the invariants II and III, and the Reynolds num- 
ber. Most of the workers have used Rotta’s hypothesis which is obtained by ignoring the 
nonlinear term. In LRR model, is taken as a constant, whereas Lumley (1978) proposed 
an expression involving II, III and the turbulence Reynolds number. This expression for /? 
was calibrated from experiments in which invariant III was negative. Gence and Mathieu 
(1980) and Choi (1983) experimentally studied the return to isotropy of homogeneous flows 
and found that the return to isotropy is slower when the invariant III is positive. Choi and 
Lumley (1984) suggested keeping the complete form of (3.27) and using their experiments 
proposed new forms of fl and 7 (see Appendix C). We will use both Lumley (1978) and 
Choi and Lumley (1984) models with the presently proposed models for the rapid term to 
see how the combinations fair. 

For the slow part of the pressure temperature-gradient correlation Shih and Lumley 
(1985)proposed 

$« = fi^ui + (3.28) 

Again most workers ignore the second term and take f\ as a constant. Shih and Lumley 
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(1986) and Shih et al (1990) use a variable form of f\ given by 


= g + (/?- 2)Ifr/6 

2 # 2 6 (//rf/3 + b tJ djj — bijdij ) 


(3.29) 


and take /2 as zero. The invariant //^ appearing in the above expression is defined by 


IU = - 4)/2 


(3.30) 


This expression is based on the joint realizability principle and the parameter H is given in 
Shih et al (1990) and f3 is the same as proposed by Lumley (1978). Both of these functions 
are given in Appendix C. In the present study, we apply the joint realizability condition 
to (3.28) keeping the full non-linear model for the pressure strain to obtain 


J €j_£ (/? + jll — 2)//rf/6 

71 2 + 02 £ (J/j/3 + M^-Mii) 

7 tl£ /J « ± Ml - Mdal + HF m 3 ' 31) 

2 (II d /3 + bijfy - b^) + d 

where /? and 7 are the same coefficients as appearing in (3.27). For 7 = 0, the Shih and 
Lumley (1986) model is recovered. In the present study we used the /? and 7 as given by 
Choi and Lumley (1984) and these are reproduced in Appendix C.The above expression 
seems to give better agreement with the experiments. In the experiments analysed, the 
value of fi varied from 6.4 - 8.0. The values used by Zeman and Lumley (1976), f\ = 7.5, 
and Newman et al (1981), f\ = 6.6, fall within this range. 


3.5 Other Models 

Other researchers, who employ the realizability concept, approach the models for 
pressure correlation terms in a similar but slightly different way. Craft et al. (1989) have 
proposed non-linear models for these correlations. They start with the similar form of 
the models as Shih and Lumley (1985), but apply realizability based only on u a u a — ► 0. 
Their models do not satisfy joint realizability and, therefore, are not fully realizable when 
a scalar is present. If we impose these conditions on their models, their coefficients must be 
changed such that their models will reduce to the same form as proposed here. Reynolds 
(1987) obtained a set of model constraints using 2D-2C conditions on X p j q i which also 
ensure realizability. His model keeps fourth powers of bij and , therefore, is more general. 
However, regardless of the complexity of the model form, the comparison between his 
model and the DNS data did not show significant improvement over the other models. 
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4. Comparison between Models and Experiments 


For model evaluation we will use the homogeneous shear flow experiment of Tavoularis 
and Corrsin (1981), DNS of the same flow by Rogers et al (1986) and the buoyant plume 
experiment of Shabbir and George (1990) and Shabbir (1991). Note that for the homo- 
geneous shear flow the buoyancy effects are negligible. For the experiments the pressure 
terms were obtained as the closing terms in the balances for the Reynolds stress and heat 
flux equations. The pressure correlations obtained in this fashion will certainly have some 
error but this is the closest we can come in obtaining these from the experiment. Also a 
direct comparison between these and their models is still a better approach in establishing 
their accuracy rather than solving the differential equations and then comparing the mean 
flow results etc. with the experiments. 


4.1 Pressure Strain Correlation 

First we look at the comparison between the pressure strain models and the experi- 
ment of Tavoularis and Corrsin (1981). The velocity field of this experiment is essentially 
the same as that of Harris et (1977) and, therefore, the results presented here are applicable 
to that experiment as well. Figure (1) shows a comparison between the experiment and 
five different models. The two complete linear models are Launder et al (1975) (LRR) and 
Zeman and Lumley (1976) and these are of the exact same form except the model coeffi- 
cients. Of the two, the Zeman and Lumley model seems to give a better overall agreement 
with the data. The Shih and Lumley (1985) model uses linear form for the slow term with 
1 3 given by Lumley (1978) and is off by a large margin. This is because j3 of Lumley (1978) 
was calibrated against flows with negative III whereas for the experiment under consider- 
ation this invariant is positive. When instead the Choi and Lumley (1984) model is used 
for the slow term, all the pressure strain components are predicted within experimental 
accuracy. Craft et al (1989) also does a good job in predicting all the components. This is 
to be expected since the coefficients in their model are determined from the homogeneous 
shear flow experiment. 

Figure (2) compares the slow term models with the DNS (19S6) results of the ho- 
mogeneous shear flow. Lumley’s (1978) linear and the Choi and Lumley (1984) nonlinear 
models give the best agreement with the data. Recall that these two models gave very 
different results for the experiment of Tavoularis and Corrsin (1981). Unlike that experi- 
ment, the Reynolds number in the DNS is much smaller. Both of these models incorporate 
the similar low Reynolds number behavior and, therefore, give similar overall results. Both 
the Zeman and Lumley (1976) and LRR (1975) model are off the data. The non-linear 
model of Craft et al (1989) performs well only for the <& xx and § xy components. Figure (3) 
compares the three different models for the rapid pressure term with the DNS of Rogers et 
al (1986). Of the three Shih and Lumley (1985) non-linear model gives the best agreement. 
The LRR (1975) model performs reasonably only for the n™ 7 "'* and W r a 2 ini components. 


17 


We also note that the Craft et. al. model performs better than the LRR model. 

Figure (4) compares the total pressure-strain models for the buoyant plume experiment 
of Shabbir and George (1990). For the streamwise component the Zeman and Lumley 
(1976) model gives the best agreement. The Shih and Lumley (1985) and LRR (1975) 
models are off the experiment near the centerline of the plume but arc within experimental 
uncertainty for the rest of the flow. Craft et. al. (1989) overprcdicts while Shih-Choi- 
Lumley underpredicts the experiment. For the radial component Zeman and Lumley 
(1976) model performs the best and, except for Shih-Choi-Lumley, all the other models 
give reasonable agreement with the experiment. For the n rj component, which is the most 
important one, Shih-Choi-Lumley, LRR (1975) and Craft et. al. (19S9) are in excellent 
agreement with the experiment. Zeman and Lumley (1976) is also within experimental 
uncertainty where as Shih and Lumley (1985) overestimates the experiment. We note that 
no single model gives good agreement with all the flows. It is our feeling that the source 
of this is either the coefficients for the slow term or the rapid term ( « 2 , given by equation 
(3.20) was assigned its limiting value). Further work is uderway in order to sort this out. 


4.2 Pressure Temperature- Gradient Correlation 

Now we compare the models for the pressure-temperature gradient correlations with 
the same three flows. First we look at the experiment of Tavoularis and Corrsin (1981) 
(figure 5). For the streamwise correlation Shih and Lumley (1985) and the present model 
provide the best agreement with the data. The Craft et. al. (1989) model overpredicts the 
data whereas Zeman and Lumley (1976) underpredicts it. The linear model of Launder 
(1975) does a reasonable job. For the correlation in the cross-stream direction the present 
model gives the best result and all the other models underpredict the experiment. Since 
this experiment was used in adjusting the constants in the present model, this comparison 
does not provide a true test for this model. 

Figure (6) shows a similar comparison with the DNS of Rogers et. al. (1986) for the 
same flow. For the streamwise correlation, again Shih and Lumley and the present model 
give the best agreement with the data. The non-linear model of Craft et. al. (1989) and 
the linear model of Launder (1975) also give reasonable results. For the correlation in 
the cross-stream direction , however, none of the models predict the data, although at the 
initial stages of the flow the Launder model does a satisfactory job. However, toward the 
last part of the simulation, where the simulation is more developed, the Launder (1975) 
model underestimates the simulation by 40%. At the same location the present model 
overestimates the data by the same amount. 

The same models are compared with the buoyant plume experiment of Shabbir and 
George (1990) in figure (7). For the streamwise component we note that the present 
model gives excellent agreement with the experiment. Both Shih and Lumley (1985) and 
Zeman and Lumley (1976) models also do a reasonable job whereas the rest of the models 
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underpredict the experiment. For the radial component, II r , the present model somewhat 
underpredicts the data near the central part of the plume, but produces the correct peak 
value and gives a very satisfactory overall agreement with the experiment. The rest of 
the models are off by a large margin. Note that the prediction of the radial component is 
extremely important for the correct calculation of the mean buoyancy field. 
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APPENDIX A 


The following axe the exact expressions for the rapid and buoyancy contributions to 
the pressure-strain and pressure scalar-gradient correlations. 

Xpjqi = -{1/**) J j Jj u q( r ) u i(r')], P jdv/\r -r'| 

X ijk = ~ (1/4tt) J j Jjuk( r )&j{ r ')],jidv/\r - r'| 

Ykij = — (1/ 47r ) J J j[uj{r)6(r% ki dv/\r -r'\ 

Y ik = -{ 1/4tt) J J J[0(r)6(r')] jik dv/\r -r'\ 

The slow terms are defined by the following expressions 

-$ij e =(1 /p)psi(ui,j + u jt i) - 2i/u >) fcU J) j t + 2e6ij/3 

e/q 2 =(1 /p)psiO,i ~(v + -y)0 ik ui ik 

The total (i.e. sum of the slow, rapid and buoyancy parts) pressure strain and pressure 
scalar-gradient correlation has been denoted by II,j and XI j respectively. 
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APPENDIX B 


The values of the model coefficients P's and 7’s as calculated from the buoyant plume 
experiment of Shabbir and George (1990) 


r/z P\ fa p z p A Pi A> Pi P^ Pa Pn 

0.00 0.47 -0.12 -2.05 -0.35 0.65 -0.20 6.40 1.95 -7.05 -1.30 

0.01 0.47 -0.12 -2.04 -0.35 0.65 -0.20 6.20 1.94 -6.85 -1.30 

0.02 0.48 -0.12 -2.02 -0.36 0.64 -0.19 5.72 1.93 -6.36 -1.28 

0.03 0.49 -0.12 -1.99 -0.37 0.63 -0.18 5.16 1.90 -5.79 -1.27 

0.04 0.50 -0.12 -1.96 -0.37 0.63 -0.17 4.66 1.88 -5.29 -1.25 

0.05 0.51 -0.13 -1.93 -0.38 0.62 -0.16 4.28 1.85 -4.90 -1.24 

0.06 0.52 -0.13 -1.91 -0.39 0.61 -0.15 4.01 1.83 -4.62 -1.22 

0.07 0.53 -0.13 -1.88 -0.40 0.60 -0.14 3.83 1.81 -4.44 -1.21 

0.08 0.53 -0.13 -1.86 -0.40 0.60 -0.13 3.73 1.80 -4.33 -1.20 

0.09 0.54 -0.14 -1.84 -0.41 0.59 -0.13 3.69 1.78 -4.28 -1.19 

0.10 0.55 -0.14 -1.83 -0.41 0.59 -0.12 3.68 1.77 -4.27 -1.18 

0.11 0.55 -0.14 -1.81 -0.41 0.59 -0.11 3.70 1.76 -4.28 -1.17 

0.12 0.56 -0.14 -1.80 -0.42 0.58 -0.11 3.74 1.75 -4.32 -1.16 

0.13 0.56 -0.14 -1.79 -0.42 0.58 -0.11 3.78. 1.74 -4.36 -1.16 

0.14 0.56 -0.14 -1.79 -0.42 0.58 -0.11 3.83 1.74 -4.41 -1.16 

0.15 0.56 -0.14 -1.80 -0.42 0.58 -0.11 3.87 1.74 -4.45 -1.16 

0.16 0.56 -0.14 -1.81 -0.42 0.58 -0.11 3.91 1.75 -4.49 - 1.17 

0.17 0.55 -0.14 -1.83 -0.41 0.59-0.12 3.93 1.77 -4.51 -1.18 

0.18 0.54 -0.13 -1.85 -0.40 0.60 -0.13 3.93 1.79 -4.53 -1.19 


r/z 7l 72 73 74 7s 

0.00 0.41 -2.82 -0.36 -0.81 -3.15 
0.01 0.41 -2.75 -0.35 -0.81 -2.96 
0.02 0.42 -2.58 -0.34 -0.80 -2.51 
0.03 0.42 -2.39 -0.32 -0.79 -1.99 
0.04 0.43 -2.22 -0.30 -0.79 -1.53 
0.05 0.43 -2.10 -0.28 -0.79 -1.19 
0.06 0.44 -2.01 -0.27 -0.81 -0.95 
0.07 0.44 -1.95 -0.27 -0.83 -0.81 
0.08 0.44 -1.92 -0.28 -0.87 -0.74 
0.09 0.44 -1.91 -0.30 -0.93 -0.72 
0.10 0.44 -1.91 -0.33 -1.00 -0.73 
0.11 0.44 -1.92 -0.37 -1.09 -0.77 
0.12 0.45 -1.94 -0.41 -1.21 -0.83 
0.13 0.45 -1.96 -0.47 -1.35 -0.88 
0.14 0.44 -1.98 -0.53 -1.53 -0.93 
0.15 0.44 -2.00 -0.59 -1.75 -0.97 
0.16 0.44 -2.02 -0.64 -2.02 -0.99 
0.17 0.43 -2.03 -0.69 -2.35 -0.98 
0.18 0.42 -2.05 -0.71 -2.76 -0.95 
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APPENDIX C 


Here we list the various model coefficients for the slow terms used in the present study. 
The coefficient /? used in the Shih and Lumley (1985) model is the one suggested by Lumley 
(1978) and is 

p =2 + exp[-7.77/R) /2 } 

{72/i?; /2 + 80.1/n[l + 62.4 (—/I + 2.3///)]} (^ + 3 III + II) 


The function H appearing in (3.29) and (3.30) is taken from Shih et. al. (1990) and is 
given by 

H = 1.1 + 0.55(/? — l)tanh[4( 1)] 

Tg 


The coefficients ft and 7 appearing in equation (3.27) and suggested by 
(1984) are given by 


P 

7 


p*F 1 / 2 

+ (1 + G X 2 ) 
p'F 1 ! 2 G 
_ (1 + Gx 2 ) i 


Choi and Lumley 


where 

£ =(//// 2) 1 ' 3 , r, = (-II/3) 1 ' 2 

X =". G = - X 1 + 0.8x 6 

p- =exp[-9.29/iJ, 1/2 l{(|§ + - [296 - 16.2(.\ + 1 )<]//} 

rCj < 


For the Zeman and Lumely (1976) model following set of coefficients were used 

1 


/? = 


(-15//) 


[1 + (1 - 24 Ilf 12 } 


/] “ ( 2 + 


Note that these revised coefficients axe taken from Zeman (19S1) and are different from 
the values originally used by Zeman and Lumley (1976). It should also be pointed out that 
the expression for P is clipped at 6.0. 
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Figure 1 . Comparison between models for II , 7 and the homoge- 
neous shear flow experiment. LRR model; Zeman and 

Lumley model (1976); — •— Shih and Lumley (1985) model;—— —Craft 

et. al. (1989) model; Shih-Choi-Lumley model; • experiment 

of Tavoularis and Corrsin (1981). 




x /h 


(Figure 1 continued from las 
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Figure 2. Comparison between the models of the slow part of the 

pressure-strain and the DNS of homogeneous shear flow. LRR 

model; Zeman and Lumley model (1976); — •— Shih and Lumley 

(1985) model; Craft et. al. (1989) model; Shih-Choi-Lumley 

model; • DNS of Rogers et. al. (1986). 





Figure 3, Comparison between the models of the rapid part of the 
pressure-strain and the DNS of homogeneous shear flow. ----- LRR 

model; Shih and Lumley (1985) model;———— Craft et. al. 

(1989) model; • DNS of Rogers et. al. (1986). 
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(Figure 3 continued from last page) 
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Figure 4. Comparison between models for Ily and the buoyant 

plume experiment LRR model; Zeman and Lumley 

model (1976); Shih and Lumley (1985) model; Craft et. 

al. (1989) model; Shih-Choi-Lumley model; e experiment of 

Shabbir and George (1990). 
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Figure 5. Comparison between models for II; and the homogeneous 

shear flow experiment Launder (1975) model; Zeman and 

Lumley model (1976); — •— Shih and Lumley (1985) model; 

Craft et. al. (1989) model; present model; • experiment of 

Tavoularis and Corrsin (1985). 
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Figure 6. Comparison between models for the II,- and the DNS of 

homogeneous shear flow Launder (1975) model; Zeman and 

Lumley model (1976); — — Shih and Lumley (1985) model; 

Craft et. al. (1989) model; present model; • DNS of Rogers 

et. al. (1986). 
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Figure 7. Comparison between models for IT,- and the buoyant 

plume experiment. Launder (1975) model; -Zeman and 

Lumley model (1976); Shih and Lumley (1985) model; 

Craft et. al. (1989) model; present model; • experiment of 

Shabbir and George (1990). 


35 



NASA 

Nation*! Aeronautics and 
Space Administration 


Report Documentation Page 


2. Government Accession No. 


1. Report No. NASA TM- 10441 3 

ICOMP-91-09; CMOTT-91-03 


4. Title and Subtitle 

Advances in Modeling the Pressure Correlation Terms 
in the Second Moment Equation 


3. Recipient's Catalog No. 


5. Report Date 

June 1991 


8. Performing Organization Code 


7. Authors) 

Tsan-Hsing Shih, Aamir Shabbir, and John L. Lumley 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135- 3191 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 - 0001 


8. Performing Organization Report No. 

E -6241 

10. Work Unit No. 

505-62-21 

1 1 . Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 

14. Sponsoring Agency Code 


1 5. Supplementary Notes 

Tsan-Hsing Shih and Aamir Shabbir, Institute for Computational Mechanics in Propulsion and Center for Modeling of 
Turbulence and Transition, Lewis Research Center (work funded by Space Act Agreement C-99066-6). Space Act 
Monitor: Louis A. Povinelli, (216) 433-5818. John L. Lumley, Cornell University, Ithaca, New York 14853. 


16. Abstract 

In developing turbulence models, different authors have proposed various model constraints in an attempt to make the 
model equations more general (or universal). The most recent of these are the realizability principle (Lumley 1978, 
Schumann 1977), the linearity principle (Pope 1983), the rapid distortion theory (Reynolds 1987) and the material 
indifference principle (Speziale 1983). In this paper we will discuss several issues concerning these principles and will 
pay special attention to the realizability principle raised by Lumley (1978). Realizability (defined as the requirement of 
non-negative energy and Schwarz’ inequality between any fluctuating quantities) is the basic physical and mathematical 
principle that any modeled equation should obey. Hence, it is the most universal, important and also the minimal 
requirement for a model equation to prevent it from producing unphysical results. In this paper we will describe in detail 
the principle of realizability, derive the realizability conditions for various turbulence models, and propose the model 
forms for the pressure correlation terms in the second moment equations. Detailed comparisons of various turbulence 
models (Launder et al. 1975, Craft et al. 1989, Zeman and Lumley 1976, Shih and Lumley 1985 and one proposed here) 
with experiments and direct numerical simulations will be presented. As a special case of turbulence, we will also 
discuss the two-dimensional two-component turbulence modeling. 


1 7. Key Words (Suggested by Author(s)) 

Turbulence 

Modeling 

Pressure correlations 


1 8. Distribution Statement 

Unclassified - Unlimited 
Subject Category 34 


1 9. Security Classif . (of the report) 

Unclassified 


20. Security Classif. (of this page) 

21 . No. of pages 

22. Price* 

Unclassified 

36 

A03 


NASA FORM 1628 OCT 86 


'For sale by the National Technical Information Service, Springfield, Virginia 22161 























